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Recently realization of TeV scale inverse seesaw mechanism in supersymmetric SO(10) framework 
has led to a number of experimentally verifiable predictions including low-mass and Z' gauge 
£SJ ■ bosons and nonunitarity effects. Using nonsupersymmetric 50(10) grand unified theory, we show 

how a TeV scale inverse seesaw mechanism for neutrino masses is implemented with a low-mass Z' 
boson accessible to Large Hadron Collider. We derive renormalization group equations for fermion 
masses and mixings in the presence of the intermediate symmetries of the model and extract the 
Dirac neutrino mass matrix at the TeV scale from successful GUT-scale parameterization of fermion 
masses. We estimate leptonic nonunitarity effects measurable at neutrino factories and lepton flavor 
violating decays expected to be probed in near future. While our prediction on the nonunitarity 
matrix element ?7 MT for degenerate right-handed neutrinos is similar to the supersymmetric 5*0(10) 
' case, we find new predictions with significantly enhanced value of its phase <5 MT ~ 10~ 4 — 10 -2 when 

partial degeneracy among these neutrino masses is adequately taken into account by a constraint 
relation that emerges naturally in this approach. Other predictions on branching ratios and CP- 
violating parameters are discussed. An important distinguishing characteristic as another test of the 
best identified minimal model is that the threshold corrected two-loop prediction on proton lifetime 
with maximum value (r p ) max . ~ 10 35 yrs. is accessible to ongoing search experiments for the decay 
Q ; p — s- e + 7r°. Simple model extensions with longer proton lifetime predictions are also discussed. 

(~| PACS numbers: 14.60. St, 12.10.Dm, 11.30.Hv, 12.15.Lk 
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I.INTRODUCTION 
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Supersymmetric grand unified theories (GUTs) provide a very attractive framework for representing particles 
and forces of nature as they solve the gauge hierarchy problem, unify three forces of nature, and also explain tiny 
neutrino masses through seesaw paradigm |l| while providing possible cold dark matter candidates of the universe. 
An evidence of supersymmetry at the Large Hadron Collider (LHC) would be a land-mark discovery which would 
certainly change the future course of physics. But, in the absence of any evidence of supersymmetry so far, it is 
worth while to explore new physics prospects of nonsupersymmetric (non-SUSY) GUTs 2-4] and, particularly, 
those based upon 5O(10) which has grown in popularity as it unifies all fermions of one generation including the 
right-handed (RH) neutrino into a single spinorial representation. It provides spontaneous origins of P (= Parity) 
and CP violations 0-0] ■ Most interestingly, in addition to predictingthe right order of tiny neutrino masses through 
mechanisms called the canonical (= type-I) [|| seesaw and type-II [§] seesaw, it has high potentiality to explain all 
fermion masses [3, 11 1 including large mixings in the neutrino sector [l2| with type-II seesaw dominance 13- TBjj . 
In fact neither seesaw mechanism, nor grand unification require supersymmetry per se. Although gauge couplings 
automatically unify in the minimal supersymmetric standard model [i"tj| , and they fail to unify through the minimal 
particle content of the standard model (SM) in one-step breaking of non-SUSY SU(5) or SO(10), they do unify 
once intermediate symmetries are included to populate the grand desert in case of non- SUSY SO(10) 0,E3fi3- In 
addition, with intermediate gauge symmetries 5O(10) also predicts signals of new physics which can be probed at 
low or accelerator energies. 

A hallmark of 5O(10) grand unification is its underlying quark-lepton symmetry Q because of which the canonical 
seesaw scale is pushed closer to the GUT scale making it naturally inaccessible to direct tests by low-energy experiments 
or collider searches. The energy scale of type-II seesaw mechanism in 5O(10) is also too high for direct experimental 
tests. In contrast to these high scale seesaw mechanisms, an experimentally verifiable and attractive mechanism that 
has been recently introduced into 5O(10) [13] is the radiative seesaw [ll| where the quark-lepton unification has no 
role to play and additional suppression to light neutrino mass prediction occurs by loop mediation proportional to a 
small Higgs quartic coupling that naturally emerges from a Plank-scale induced term in the GUT Lagrangian. The 
model predicts a rich structure of prospective dark matter candidates also verifiable by ongoing search experiments. 
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It has been further noted that this embedding of the radiative seesaw in 5*0(10) may have a promising prospect 
for representing all fermion masses. A number of other interesting neutrino mass generation mechanisms including 
type-Ill seesaw, double seesaw, linear seesaw, scalar- triplet seesaw have been suggested and some of them are also 
experimentally verifiable [l[. 

In the context of non-SUSY 50(10) in this work our purpose is to explore the prospects of another neutrino mass 
generation mechanism called the inverse seesaw 22j which is also different from canonical or Type-II seesaw mechanism 
and has the potentiality to be experimentally verifiable because of the low scale at which it can operate although 
higher scale inverse seesaw models have been suggested 23, 3|- I n a large class of models 2^- 28|, the implementation 
requires the introduction of fermionic singlets under the gauge group of the model. Likewise, its implementation in 
50(10) introduces a new mass scale /xg into the Lagrangian corresponding to the mass matrix of the additional singlet 
fermions of three generations and the TeV-scale seesaw requires this parameter to be small. There is an interesting 
naturalness argument in favor of its smallness based upon exact lepton number conservation symmetry 25, |29j . 
Below the TeV scale, in the limit (is 0, the corresponding Lagrangian has a leptonic global U(l) symmetry which 
guarantees left-handed neutrinos to remain massless. The small value of /is essentially needed to match the neutrino 
oscillation data with the TeV-scale inverse seesaw formula may be taken as a consequence of very mild breaking of 
the exact global symmetry. Thus the small value of the parameter protected by the exact lepton number conservation 
is natural in the 't Hooft sense (30j . However, inspite of such interesting naturalness argument, there has been no 
dynamical understanding of its origin so far although an interpretation using Higgs mechanism has been given in 
the context of a model with extended gauge, fermion and Higgs sectors 31 1. In a different class of 5O(10) models, 
the small singlet fermion mass parameter has been generated radiatively 32j where more nonstandard fermions have 
been found to be necessary. In SUSY 5O(10) singlet fermions have also been used to derive new forms of fermion 
mass matrix while predicting standard fermion mass ratios [33^ and to obtain new seesaw formula for neutrinos while 
explaining baryon asymmetry of the universe through leptogenesis [34| . While most of the inverse seesaw models 
need gauge singlet fermions under the SM gauge group or its extensions 22-2^, 31, 32| and the use of 5O(10)-singlet 
fermions may point to the disadvantages of the corresponding GUT-based models, extended electroweak theory based 
upon SU(3)l x SU(3)n, x U(l) gauge symmetry [35j contains such singlets in its fundamental representations. To give 
some examples of GUTs, one 5O(10)-singlet fermion per generation is automatically contained in the 27— dimensional 
fermion representation of E 6 ^ where 27 = 16 + 10 + 1 under 5O(10) but 27 = (3, 3*, 1) + (3, 1, 3*) + (1, 3, 3*) under 
SU(3)l x SU(3)b, x SU(3)c and, in the latter case, an additional discrete Z 3 symmetry is needed to qualify it as 
a trinification GUT model (37J. Interesting properties of 5C/(3) 3 gauge theory including experimentally verifiable 
predictions at accelerator energies have been discussed [HI, 39] . Gauge boson mediated proton decays are suppressed 
in 5C/(3) 3 type of models. In addition to the RH neutrino and the other singlet fermion needed for inverse seesaw, 
these models ( SU(3) 3 and Eq) also contain 10 nonstandard fermions per generation and no experimental data are yet 
available on their masses at low energies so as to persue the present bottom-up approach to derive the Dirac neutrino 
mass matrix from fermion mass fits at the GUT scale. The same argument holds against any other model that may 
contain additional nonstandard fermions beyond the RH neutrino and the singlet-fermion needed for inverse seesaw. 

Regarding the potentiality of SO(10) motivated inverse seesaw in the visible sector, the same quark-lepton symmetry 
that forces canonical seesaw scale to be far beyond the experimentally accessible range, makes the TeV-scale inverse 
seesaw predict observable nonunitarity effects as new physics signals verifiable at low and accelerator energies and at 
neutrino factories (29| . 

Recently in a series of interesting investigations, using inverse seesaw mechanism, Bhupal Dev and Mohapatra [29j] 
have shown that SUSY 5O(10), besides admitting a low spontaneous breaking scale of SU(2)l x SU(2)r x U(1)b-l x 
SU(3)c(= G2213) gauge symmetry with right handed gauge bosons and Z' accessible to LHC, is also capable of 
fitting all fermion masses and mixings at the GUT scale while predicting observable nonunitarity effects. The model 
has been also shown to account for the observed baryon asymmetry of the universe through leptogenesis caused due 
to the decay of TeV scale masses of the pseudo Dirac RH neutrinos [I(| ■ Currently considerable attention has been 
devoted to propose models with an extra neutral Z' gauge boson which may also emerge from Pati-Salam or left-right 
gauge theories, or 5O(10) and Eq grand unified theories, or also from string inspired models 41] . 

Different from the supersymmetric 5O(10) model of ref.[29], here we adopt the view that there may not be any 
manifestation of supersymmetry at accelerator energies and that the actual parity restoration scale may be high. 
Instead of both the and the Z' boson masses being low, there may be only some remnants of high scale left-right 
symmetry or quark lepton symmetry manifesting at low and accelerator energies as smoking gun signatures such 



as the Z' [4JJ, |42[ gauge boson and the associated nonunitarity effects of the TeV-scale inverse seesaw. With this 
point of view in this work we show that a non-SUSY 5O(10) with SU(2)l x U(1)r x U(1)(b-l) x SU(3)c(= G2113) 
gauge symmetry at the TeV scale and left-right gauge theory at higher intermediate scale, with or without D-parity, 
achieves precision gauge coupling unification, and predicts a low mass Z' making them suitable for implementation of 
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TeV-scale inverse seesaw mechanism. The model can also be verified or falsified through its predictions on observable 
nonunitarity effects and additional contributions to lepton flavor violations. Another testing ground for the model 
could be through the SO(10) prediction on gauge boson mediated proton decay on which dedicated search experiments 
are ongoing. 

We derive renormalization group equations in the presence of two intermediate gauge symmetries for running 
fermion masses and mixings, and determine the Dirac neutrino mass at the TeV scale from successful fits to the 
fermion masses at the GUT scale. In this approach we find a simple relation between the RH neutrino masses in the 
model. We also point out a different type of relation in the partial degenerate case that permits much lower values of 
RH neutrino masses resulting in a CP-violating phase increased by 2 — 4 orders larger than the degenerate case. Some 
of our predictions include branching ratios for /i — > ej enhanced by 1 — 2 orders. Out of the two minimal models, 
while the intermediate scale D-Parity conserving model is ruled out by proton decay constraint, the proton lifetime 
for p — >• e + 7r° in the intermediate scale D-Parity nonconserving model is predicted to be well within the accessible 
range of ongoing search experiments. We have also discussed simple extensions of the two models with longer proton 
lifetime predictions. This method can also be implemented using Pati-Salam model or left-right models @Jl]. 

This paper is organized in the following manner. In Sec. 2 we briefly discuss the model and carry out gauge 
coupling unification and proton lifetime predictions in Sec. 3. With a brief explanation of inverse seesaw mechanism in 
Sec. 4 we summarize relevant formulas encoding nonunitarity effects and lepton flavor violations. In Sec. 5 we discuss 
renormalization group evolution of fermion masses and mixings to the GUT scale in the presence of nonsupersymmetric 
gauge theories G2113 and G2213 at intermediate scales. In this section we also show how fermion masses are fitted 
at the GUT scale and information on the Dirac neutrino mass matrix is obtained. Nonunitarity effects are discussed 
in Sec. 6 with predictions on the moduli of relevant matrix elements. In Sec. 7 we give predictions on CP-violating 
parameters and lepton flavor violation where we also discuss possible limitations of the present models. In Sec. 8 
we provide a brief summary and discussion along with conclusion. In the Appendix A we provide beta function 
coefficients for gauge coupling unification while in Appendix B we summarize derivations of renormalization group 
equations (RGEs) for fermion masses and mixings. 



II.THE MODEL 
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There has been extensive investigation on physically appealing intermediate scale models ll 
SUSY 50(10). Although in the minimal two step-breaking of non-SUSY 50(10) models [2 
chain with a sufficiently low scale to implement the inverse seesaw, the following chain with two intermediate gauge 
symmetries appears to be quite suitable, 
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SU{2) L x U{1) R x U(1) B -l x SU(3) C [G 2113 ] 

SU(2) L x U(1) Y x SU{3) C [SM] 

SU(3) C x U(1) Q , 
where we will consider two possibilities for G/. 



(1) 



As Model-I G/ = SU(2) L x SU(2) R x U(1) B -l x SU(3) c [= G2213K32L ^ 92b) is realized by breaking the GUT- 
symmetry and by assigning vacuum expectation value (VEV) to the D-Parity odd singlet in 45// @ • As the left-right 
discrete symmetry is spontaneously broken at the GUT scale, the Higgs sector becomes asymmetric below fj, = Mjj 
causing inequality between the gauge couplings giL and g2R ■ This model does not have the cosmological domain wall 
problem. The second step of symmetry breaking takes place by the right-handed (RH) Higgs triplet <tr(1, 3, 0, 1) C 45// 
whereas the third step of breaking to SM takes place by the G2ii3-submultiplet x5f (l, 1/2, —1/2, 1) contained in the 
RH doublet of 16//. It is well known that SM breaks to low energy symmetry by the SM Higgs doublet contained 
in the bidoublet (2,2,0,1) under G2213 which originates from 10/r of 5O(10). This is the minimal particle content 
for the model to carry out the spontaneous braking of GUT symmetry to low-energy theory. But a major objective 
of the present work is to explore the possibility of observable nonunitarity effects for which it is required to extract 
information on the Dirac neutrino mass matrix (Mb) from a fit to the fermion masses at the GUT scale and this 
is possible by including two Higgs doublets instead of one [2^. We assume these doublets to originate from two 
separate bidoublets contained in 10^ (a = 1,2). Implementation of inverse seesaw also requires the minimal extension 
by adding three 5O(10)-singlet fermions Si(i = 1,2,3), one for each generation [221 ]. 
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As Model-II, we treat the GUT symmetry to be broken by the VEV of the G 2 2i3-singlet (1, 1, 0, 1) C 210 H which 
is even under D- parity @. This causes the Higgs sector below GUT-scale to be left-right symmetric resulting in 
equal gauge couplings in Gj = G2213d{92L — 92R)- For the sake of simplicity we treat the rest of the symmetry 
breaking patterns of Model-II similar to Model-I and we assume the presence of three singlet fermions. We call these 
two models, Modcl-I and Model-II, as minimal models with two low scale Higgs doublets in each. We now examine 
precision gauge coupling unification for these two models. 



III.GAUGE COUPLING UNIFICATION AND PROTON LIFETIME 



In this section we examine gauge coupling unification in the minimal Model-I and and the minimal Model-II and 
make predictions on proton lifetimes while we also predict the corresponding quantities in their simple extensions. 

A. Unification in minimal models 

It was shown in that with G2113 gauge symmetry at the lowest intermediate scale in 5*0(10) there is substantial 
impact of two-loop effects on mass scale predictions in a number of cases. The one-loop and the two-loop beta-function 
coefficients for the evolution of gauge couplings 44, 4|| for Model-I and Model-II with two Higgs doublets for each 
case are given in Appendix A. We have also included small mixing effects [43 , 46 1 due to two abelian gauge factors 
U(l) R xU(l) (B -L) in both the models below the M+ scale. Using sin 2 9 W (M Z ) = 0.23116±0.00013, a~ 1 {M z ) = 127.9 
and as(Mz) = 0.1184 ± 0.0007 [47j we find that with Mz> ~ M R o ~ 1 TeV precision unification of gauge couplings 
occurs for the following values of masses at one-loop and two-loop levels for the Model-I, 



Mu = 
Mu = 



10 i 5 .978 GeViM oi + = 10 io.787 GeV; Q oj = q.02253, (one - loop), 
10 i5.530 GeVj Mr+ = 10 n.i5 GeVj ac = _o2290 (two - loop). 



(2) 



The RG evolution of gauge couplings at two-loop level is shown in Fig. [T] exhibiting precision unification at Mjj = 
10 15 53 GeV. In Model-II coupling unification occurs with similar precision but at Mjj = 10 1517 GeV. 
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FIG. 1: Gauge coupling unification in Model-I with two-loop values Mu = 10 15 53 GeV and M R+ = 10 1115 GeV with a low mass Z 1 

gauge boson at M R o ~ 1 TeV. 



The decay width of the proton for p — > e + n is [48 1 
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■^qIS.978 


^q10.787 


^glS.530 


■^gll.150 


43.67 


2.726 


1.08 x io 36±0 - 32 


2 x io 34±0 - 32 


II 


^q15.56±0.08 


10 11.475 


^q15.17±0.08 


10 11.750 


42.738 


2.670 


2.44 x io 34±0 - 32 


6.3 x io 32±0 ' 32 



TABLE I: GUT scale, intermediate scale and proton lifetime predictions for nonsupersymmetric SO(10) models with TeV scale Z' boson 
and two Higgs doublets as described in the text. The uncertainty in the proton lifetime has been estimated using 3<r uncertainty in 



where R = [{A 2 SR + + |Kd| 2 ) 2 ] for 50(10), V ud = 0.974 = the (1,1) element of Vckm for quark mixings, 

Asl{Asr) is the short-distance renormalization factor in the left (right) sectors and Al = 1.25 = long distance 
renormalization factor. Mjj — degenerate mass of 24 superheavy gauge bosons in 5*0(10), an = hadronic matrix 
element, rn p = proton mass = 938.3 MeV, — pion decay constant =139 MeV, and the chiral Lagrangian parameters 
are D = 0.81 and F = 0.47. With a H = a H (l + D + F) = 0.012 GeV 3 estimated from lattice gauge theory 
computations, we obtain An ~ AlAsl — AlA$r — 2.726 for Model-I. The expression for the inverse decay rates for 
both the minimal models is expressed as, 

r-i(p-e^) = (1.01 x lQ3*Yrs) ( °-° 12 ^V ( gjggV ( iZ^V 



^F q J V2.98 x 10 15 GeV, 

where the factor F q — 2(1 + |Vud| 2 ) 2 — 7.6 for S'O(IO). Now using the estimated values of the model parameters in 
each case the predictions on proton lifetimes for both models are given in Table U where the uncertainties in unification 
scale and proton lifetime have been estimated by enhancing the error in to 3tr level. It is clear that with maximal 
value {T p ) max . = 7 x 10 34 Yrs., Model-I predicts the proton lifetime closer to the current experimental lower bound 
{T p ) expt (p — > e + 7r°) > 1.2 x 10 34 Yrs. (401 which is accessible to ongoing proton decay searches in near future [H(|. On 
the other hand Model-II is ruled out at two-loop level as it predicts lifetime nearly two orders smaller. The reduction 
of lifetime by nearly two-orders compared to one-loop predictions in both cases is due to the corresponding reduction 
in the unification scale by a factor of ~ 1/3. 

The fact that the Model-I admits a low (B — L) breaking scale corresponding to a light Z' accessible to accelerator 
searches makes this non-SUSY model suitable to accommodate inverse seesaw mechanism. Unlike the SUSY 5O(10) 
model [lij], here the W R bosons are far beyond the LHC accessible range. 

B. Unification in simple model extensions 

Although the minimal Model-I clearly satisfies the proton decay constraint to accommodate TeV scale seesaw, 
we study simple extensions of both models to show that they can evade proton lifetime constraint in case future 
experiments show t p to be substantially longer than 10 35 Yrs. We use an additional real color octet scalar Cs(l, 0, 8) C 
45# where the quantum numbers are under the SM gauge group and allow its mass to vary between 1 TeV and the 
GUT scale. Making it light would require additional fine tuning of parameters. Recently such a light scalar has been 
used in models with interesting phenomenological consequences and if the particle mass is in the accessible range, it 
may be produced at LHC with new physics signatures beyond the standard model [HlT ]. 

The presence of this scalar octet with lower mass makes the evolution of a^, 1 (jjl) flatter thereby pushing the GUT 
scale to higher values. In Fig [5] we plot predicted proton lifetimes in the extended G2213 and G2213D models as a 
function of the octet mass m%. It is clear that such a simple extension of the two models can easily satisfy proton 
lifetime requirements in future experimental measurements even if they are found to be much longer than the current 
limit. 

Then while the minimal Model-I can be easily chosen for inverse seesaw, both the models with such simple extension 
and possessing TeV scale U(1)(b-l) breaking scale qualify for the same purpose. 



IV.INVERSE SEESAW AND FORMULAS FOR CP AND LEPTON FLAVOR VIOLATIONS 



For the phenomenological study of nonunitarity effects we confine to the Model-I and all our analyses are similar 
for Model-II. Introducing additional 5O(10)-singlet fermions (S) for three generations, the Yukawa Lagrangian at the 
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FIG. 2: Variation of proton lifetime as a function of color octet mass in simple extensions of Model-I (double dot-dashed line) and 
Model-II (dashed line). The horizontal solid line with error band is the prediction of the minimal Model-I while the horizontal 

dot-dashed line is the experimental lower bound for p — > e + ir°. 



GUT scale gives rise to the effective Lagrangian near the second intermediate scale /i = M R ~ l TeV, 

£ Yuk = Y a 16 ■ 16 • 10 H a + y x 16 • 1 • 16 H f + Msl ■ 1 

D (Y a ^ L iJj R <l> a + y x ^ R Sx° R + H.c.) + S T fi s S, 



(5) 



where the first (second) equation is invariant under 50(10) (G2113) gauge symmetry. The left-handed (LH) and 
the right-handed (RH) fermion fields ^,(2,0,-1/2,1), (1, 1/2, — 1/2, 1) with their respective quantum numbers 
under G2113 are contained in the spinorial representation 16 C SO(10) and the two Higgs doublets $(2, ±1/2, 0, 1) C 
10 h C 5O(10). The Lagrangian has a new mass scale /zg corresponding to the mass matrix of the 5O(10)-singlet 
fermions. Denoting the RH neutrino mass as Mr = y x v x where v x — (Xr) and the Dirac mass matrix for neutrino 
as Md = Y v v u where v u is the VEV of the up-type Higgs doublet, eq.JS]) gives the mass part of the neutrino sector 
in the Lagrangian in the flavor basis after the symmetry breaking G2113 — > SM 



Anass = (DM D N + NM R S + H.c.)+S T fi s S, 



which, in the (y, N, S)l basis, leads to a mass matrix [22J, |2 




(6) 



(7) 



Denoting X = M^M^ 1 , block diagonalization of eq.© under the condition Mr 3> Mjj 3> jjig leads to the inverse 
seesaw formula for light neutrino mass matrix, 



MoM-^siMlyHll = X[i s X 7 



(8) 



It is clear that the TeV-scale inverse seesaw formula is tenable and appropriate to fit the light neutrino masses 
provided /is is the smallest of the three mass scales occurring in eq.([7]). Based upon symmetry, there exist interesting 
naturalness arguments in the literature in favor of smallness of /is. In the limit /ig -> a leptonic U(l) global 
symmetry is restored in the Lag rangian signifying exact conservation of lepton number that guarantees left-handed 
neutrinos to be massless [25, 29[ 34 1. In particular, a small and nonvanishing value of (is can be viewed as a slight 



breaking of the global U(l) symmetry. Thus the smallness of /is, desired in the TeV-scale inverse seesaw mechanism, 
which is protected by the global symmetry in the 't Hooft sense [3pj ] . is natural even though there is no dynamical 
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understanding for such a small parameter. This view for the naturally small parameter fi$ being followed in the 
present work has been adopted in [29j j and by a number authors earlier pursuing inverse seesaw mechanism 25 1 
although its interpretation through Higgs mechanism has been discussed in a model with extended gauge, fermion 
and Higgs sectors [31J and possibility of its radiative origin has been explored 32j. 

The physics underlying nonunitarity effects have been discussed at length in several recent papers 52-59J where 
relevant formulas have been utilized. Although the PMNS matrix U diagonalizes the light neutrino mass matrix of 
three generations where 

U^myU* = diag(mi,TO2,m 3 ) = rh u , (9) 

the appropriate diagonalizing mixing matrix for the inverse seesaw matrix of eq.® is a 9 x 9 matrix V, 

V^MuV* = M = diag(m i ,rn Rj ,m ilk ), 

{i,j,k = 1,2,3), (10) 



and this can be expressed in block partitions, 



v = r 3x3 JH = (ii) 

where the nonunitary Vax3 matrix now represents the equivalent of the full PMNS matrix, 

N = V 3x3 ~ (1 - ~XX*)U 

~ (l- V )U. (12) 

Denoting the corresponding nine component eigenstate as (z)j, iVj, Ni) T , the six component heavy eigenstate as P T = 
(N 1 ,N2,N 3 ,Ni,N 2 ,N 3 ) T and JC = V 3X 6 — (0,X)y 6x6 , in the leading order approximation in X, the light neutrino 
flavor eigenstate and the charged current Lagrangian in the mass basis are, 



v 



T = Mv T + JCP 1 



C CC = -^j=l L -fvW- + H.c. 

~ - 9 -lLl Ll »(Mv T + JCP T )W- +H.c. (13) 

The parameter rj = XX^ /2 characterizing nonunitarity of the neutrino mixing matrix can have dramatic impact on 
leptonic CP-violation and branching ratios for processes with lepton flavor violation (LFV), 

c J + AJ«, (14) 
where J is the well known CP- violating parameter due to unitary PMNS matrix U 

J = cos #i2C0S 2 #i3 cos #23 sin #12 sin #i3 sin #23 sin <5, (15) 
and the nonunitarity contributions are, 

+ VPfUaiU-fjUajUpt 

+ TlfrUaiUpjU^U*;). (16) 

Very recently sin #13 has been measured 60] to be small and nonvanishing although no experimental information is 
available on the leptonic CP-phase S. Even in the limiting case of vanishing unitarity CP-violation corresponding to 
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sin 6*i3 —> 0, or S —> 0, tt for nonvanishing #13 nonunitarity effects caused due to r\ may not vanish. In the modified 
charged current interaction in eq. (|13p . the heavy neutrinos contribute to LFV decays with branching ratios 61 j 



BR{l a -> l pl ) 



256vr 2 M4r Q 



6 



in 



Ri 



I{x) 



2x 3 + 5x 2 



x 



3x In x 



4(1 -a;) 3 2(1 - x) 4 



(17) 



In ea.(|17p the total decay width r a for lepton species l a with lifetime r Q is evaluated using r a = — where r M = 
(2.197019 ± 0.000021) x 10~ 6 sec and t t = (290.6 ± 1.0) x 10~ 15 sec. 

The matrix element {K,lO) a p oc rj a p may lead to significant lepton flavor violating (LFV) decays in the TeV scale 
seesaw whereas LFV decays are drastically suppressed in Type-I seesaw in SO(10). The procedure for estimating 
these effects has been outlined in [2!| which we follow. The Dirac neutrino mass matrix at the TeV scale which we 
derive in the next section is central to the determination of nonunitarity effects. 



V.RG EVOLUTION OF FERMION MASSES AND DETERMINATION OF Mn 

The determination of the Dirac neutrino mass matrix Mp(Mro) at the TeV seesaw scale is done in three steps 
[29| : (1) Derivation of RGEs for the specific model and extrapolation of masses to the GUT-scale, (2) Fitting the 
masses at the GUT scale and determination of Md(Mcut), (3) Determination of Md(Mro) by top-down approach. 



A. RGEs and extrapolation to GUT scale 

At first RGEs for Yukawa coupling matrices and fermion mass matrices are set up from which RGEs for mass eigen 
values and CKM mixings are derived in the presence of G2113 and G2213 symmetries. RGEs in dynamical left-right 
breaking model has been derived earlier (62[. 

Denoting <&i j2 as the corresponding bidoublets under G2213 they acquire VEVs 

< $1 > = 

< $2 > - [ " " I • (18) 
Defining the mass matrices 

M u = Y u v u , M D = Y u v u) M d = Y d v d , 

M e = Y e v d ,M R = y x v x , (19) 

we have derived the new RGEs in the presence of non-SUSY G2113 and G2213 gauge symmetries for matrices Y it M i7 i — 
u, d, e, N, the mass eigenvalues TOj, i = u,c, t, d, s, b, e, /1, r, N±, N%, N3, and the CKM mixing matrix elements as given 
in the Appendix B. We use the input values of running masses and quark mixings at the electroweak scale as in 
refs. [tEEI and the resulting CKM matrix with the CKM Dirac phase S q = 1.20 ± 0.08 




V, 



CKM 



/ 0.9742 0.2256 0.0013-0.0033^ 

-0.2255 + 0.000U 0.9734 0.04155 | . (20) 

^ 0.0081 - 0.0032i -0.0407 - 0.0007i 0.9991 



We use RGEs of the standard model for /j, = Mz to M R — 1 TeV. With two Higgs doublets at /1 > M R we use the 
starting value of tan^ = v u /v d = 10 at n = 1 TeV which evolves to reach the value tan/3 ~ 6.9 at the GUT scale. 



Using the bottom- up approach discussed earlier 63j and the RGEs of Appendix B, the resulting quantities including 
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the mass eigen values mi and the Vqkm at the GUT scale are 64 1 
H = M GUT : 



m e = 0.48 McV,m M = 97.47 MeV,m T = 1.8814 GeV, 
m d = 1.9 MeV,m s = 38.9 MeV,m fc = 1.4398 GeV, 
m u = 1.2 MeV, m c = 0.264 GeV, m t = 83.04 GeV, 



(21) 



0.9748 0.2229 -0.0003 - 0.0034i N 

Vcku{M gut ) = | -0.2227 - O.OOOli 0.9742 0.0364 . (22) 

0.0084 - 0.0033i -0.0354 + 0.0008i 0.9993 j 

B. Determination of Md 

With Higgs representations 45#, 16# , 10jj, the dim. 6 operator [2§| 

^16 i 16jl0 H 45 H 45H, (23) 

with M ~ Mpi or M ~ M string , is suppressed by (Mu/M) 2 ~ 10~ 3 - 10~ 5 for GUT-scale VEV of 45 H and acts as 
an effective 126 h operator to fit the fermion masses at the GUT scale where the formulas for mass matrices are 

M u = G u + F, M d = G d + F, 

M e = G d — 3F, Md — G u —3F. (24) 

In ea.(|24|) the matrices Gk — Yfc.16.16 < 10% >,k — u,d and F are derived from eq. ([23|) . Using a charged-lepton 
diagonal mass basis and eq. (|2"Tj) and eq. (|24| we get, 

M e (M GUT ) = diag(0.0005, 0.098, 1.956) GeV, 

Gd.ij = SFij, (25) 

Assuming for the sake of simplicity that the matrix F is diagonal leads to the conclusion that the matrix Gd is also 
diagonal. This gives relations between the diagonal elements which, in turn, determine the diagonal matrices F and 
Gd completely 

Gd,ii + Fu = mi, (i = d,s,b), 
Gd,jj~^Fjj = mj, (j' = e,/i,r), (26) 

F = diag- (ma - m e ,m s - m^mb — m T ), 

= diag(3.75 x 10~ 4 , -0.0145, -0.3797) GeV, 
G d = diagi(3?rid + m e ,3m s + m ll ,'&mb + m T ), 

= diag(0.0016, 0.0544, 1.6709) GeV, (27) 

where we have used the RG extrapolated values of eq. (l2TT) . Then using eq. (l24l) . eq. (l27l) and the assumed basis gives 
the mass matrices M u and G u , 



M U (M GUT ) 



( 0.0153 0.0615- 0.0112i 0.1028- 0.2706^ 
0.0615 + 0.0112* 0.3933 3.4270 + 0.0002i ] GeV, (2M 
y0.1028 + 0.2706i 3.4270- 0.0002 82.90 
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' 0.0150 0.0615- 0.0112i 0.1028 - 0.2706i\ 
G U (M GUT ) = 0.0615 + 0.0112i 0.4079 3.4270 + 0.0002z GeV. 

v 0.1028 + 0.2706i 3.4270 - 0.00022' 83.01 / 

Now using eq. (f27f and eq. ([29|) in eq.(|24|) gives the Dirac neutrino mass matrix Mb at the GUT scale 

/ 0.0139 0.0615- 0.0112* 0.1029- 0.2707^ 
M D (M GUT ) = 0.0615 + 0.0112i 0.4519 3.4280 + 0.0002z' I GeY. 
^0.1029 + 0.2707i' 3.4280- 0.00022 83.340 



(29) 



(30) 



We then use the RGE for Mb given in Appendix A to evolve Mb{Mqut) to Mb(M r +) and then from Md(M r +) to 
Md(M r o) in two steps and obtain, 



M D (M R o) 



( 0.0151 0.0674- 0.0113? 0.1030 - 0.27182 > 
0.0674 + 0.0113i 0.4758 3.4410 + 0.0002*' | GcY. 
y0.1030 + 0.2718«' 3.4410 - 0.00022' 83.450 



(31) 



VI.NONUNITARITY DEVIATIONS IN LEPTON MIXING MATRIX 



From eq. (]T2j) it is clear that any nonvanishing value of rj is a measure of deviation from the unitarity of the PMNS 
matrix. Using the TeV scale mass matrix for Mb from eq. (|31[> and assuming 



M R = diag(m Rl ,m R2 ,m R3 ), 



(32) 



results in 



V 



\x.X^ = M D M R 2 M D /2, 



M Dak M* T 



(33) 



fe=l,2,3 H k 

For the sake of simplicity assuming degeneracy of RH neutrinos masses m R = m Ri (i = 1, 2, 3) gives, 

/ 0.0447 0.1937-0.4704* 4.4140-11.360^ 
0.1937 + 0.4704i 6.036 144.40 - 0.00022' 

^4.4140+11.3602 144.40 + 0.0002z' 3488.0 



lGcV 2 



n 



(34) 



The deviations from unitarity in the leptonic mixing is constrained, for example, by deviations from universality 
tests in weak interactions, rare leptonic decays, invisible width of Z boson and neutrino oscillation data. The bounds 
derived at 90% confidence level from the current data on the elements of the symmetric matrix are summarized in 



\Vrr\ < 2.7 x HT 3 , |?7 W | <8.0x HT 4 , 
\ Vee \ < 2.0 x lO- 3 , < 3.5 x 10" 5 , 

|ry eT | < 8.0 x 10~ 3 , |^ ir | < 5.1 x 10~ 3 . 



(35) 



In the degenerate case the largest element in eq. (|34[) when compared with \rj TT \ of eq. (|35[) gives the lower bound on 
the RH neutrino mass, 



m R > 1.1366 TeV, 



(36) 



which is only 7% higher than the SUSY 5O(10) bound (m R )susY > 1-06 TeV [29[. Using this lower bound for other 
elements in ea. (l3"4l yields 



\Vw\ < 4.672 x lO" 6 , 
\rj ee \ < 3.460 x 10~ 8 , < 3.938 x 10 

\r) eT \ < 9.436 x 10~ 6 , \r)„ T \ < 1.1178 x 10~ 4 . 
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(37) 
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m Rl2 
(GeV) 


m,R 3 
(GeV) 


m Ri,2 
(GeV) 


m,R 3 
(GeV) 


48.0 


5572.83 


500.0 


1140.66 


50.0 


3324.69 


600.0 


1139.11 


100.0 


1286.51 


700.0 


1138.18 


150.0 


1195.81 


800.0 


1137.57 


200.0 


1168.32 


900.0 


1137.16 


300.0 


1149.80 


1000.0 


1136.87 


400.0 


1143.53 


1136.58 


1136.58 



TABLE II: Variation of third generation RH neutrino mass mji a as a function of first or second generation RH neutrino mass in the 
partially degenerate case = nin 2 predicted by nonunitarity through nonsupersymmetric 50(10). 



As in SUSY SO(10) [29[ , these predicted bounds are several orders lower than the current experimental bounds and 
they might be reached provided corresponding LFV decays are probed with much higher precision. But compared to 
SUSY 5O(10), in this model the upper bound is nearly 2 times larger for |ry MT |, 3 times larger for |r/ MAt |, and nearly 
40% smaller in the case of \r) eT \. It is interesting to note that in the the present non-SUSY 5O(10) model while some 
of the nonunitarity effects are comparable to the results of [2i| , others are distinctly different as shown in the next 
section. 

We note in this model that when RH neutrino masses are nondegenerate, they are also constrained by the experi- 
mental lower bound on rj TT and the corresponding relation obtained by saturating the bound is 



0.0845 



11.8405 + 6963.9 

m R 2 m R 3 



= 2.7 x 10" 



(38) 



where the numerators inside the square bracket are in GeV 2 . Using partial degeneracy, m Rl = tor 2 ^ TnR 3 leads to 
the relation between the RH neutrino masses as given in Table [TT1 A plot of mn 3 vs. m Ri (i = 1, 2) is shown in Figj3] 
exhibiting increase of mn 3 with decrease of m Ri . The two asymptotes in the hyperbolic curve are at = mn 2 ~ 47 
GeV and m R , ~ 1136.6 GeV. 



> 
a 
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2000 



1000, 




100 150 200 
m R (GeV) 



300 



FIG. 3: Variation of the third generation RH neutrino mass mn 3 as a function of first or second generation neutrino mass mn 1 or mn 2 

in the partially degenerate case for which = mn 2 . 



VII.ESTIMATIONS OF CP AND LEPTON FLAVOR VIOLATIONS AND DISCUSSIONS 

Two important physical applications of inverse seesaw are leptonic CP and flavor violation effects reflected through 
the elements, both moduli and phases, of the ?/-matrix and the relevant formulas have been discussed in Sec. 4. The 
inverse seesaw formula of eq.® has three matrices out of which Mn has been determined by fitting the charged 
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mR t — m,R 2 
(GeV) 


niR 3 
(GeV) 




Sen 


\Ver\ 


Set 


I??Mt| 




1136 


1136 


3.938 x 10" 7 


1.180 


9.436 x 10 -6 


1.20 


1.118 x 10~ 4 


1.3 x 10" 6 


500 


1141 


4.222 x 10~ 7 


1.071 


9.576 x 10" 6 


1.166 


1.136 x 10~ 4 


2.0 x 10" 4 


100 


1286 


1.848 x 10~ 6 


0.308 


1.687 x 10~ 5 


0.563 


1.691 x 10~ 4 


5.0 x 10~ 3 


50 


3325 


6.733 x 10 -6 


0.172 


4.806 x 10~ 5 


0.202 


3.424 x 10~ 4 


1.0 x 10~ 2 



TABLE III: Predictions of moduli and phases of nonunitarity parameters as a function of RH neutrino masses. 



m -Ri,-R 2 

(GeV) 


m R3 
(GeV) 


BR(n ->• e 7 ) 


BR(t -> &y) 


BR(t fij) 




AJ 23 


AJ, 23 


AJ- M 31 


ax 1 ! 


1136 


1136 


2.0 x 10" 16 


2.1 x 10" 14 


3.0 x 10~ 12 


-1.3 x 10" 6 


-1.6 x 10~ 6 


1.6 x 10" 6 


1.6 x 10" 6 


4.0 x 10" 6 


500 


1140 


2.0 x 10 -16 


1.9 x 10" 14 


2.7 x 10~ 12 


-1.3 x 10 -6 


-1.6 x 10 -6 


1.6 x 10 -6 


1.6 x 10 -6 


4.0 x 10 -6 


100 


1286 


1.4 x 10" 15 


2.2 x 10~ 14 


2.2 x 10 -12 


-1.2 x 10~ 6 


-1.6 x 10~ 6 


2.2 x 10" 6 


1.3 x 10~ 6 


4.1 x 10" 6 


50 


3325 


1.1 x 10" 14 


1.1 x 10~ 13 


5.5 x 10~ 12 


-1.0 x 10~ 6 


-1.8 x 10 -6 


4.1 x 10~ 6 


7.4 x 10 -7 


4.3 x 10 -6 



TABLE IV: Nonunitarity predictions of leptonic CP-violating parameters and branching ratios for lepton flavor violating decays fj, — > ej, 

t — > e7, and r — > fi-y as a function RH neutrino masses 



fermion masses and mixings, but since the other two matrices, Mjv and mus, can not be completely determined by 
using the neutrino oscillation data alone, we make plausible assumptions. In addition to the fully degenerate case we 
also examine consequences of partial degeneracy with = mn 2 . 

From eq.([8]), the nommitary PMNS matrix J\f = (1 — rfjU and the relation m v = Mfh v M T give, 

lis = X- l Nm u N T {X T )- 1 . (39) 

We construct the unitary matrix U using standard parameterization, 



C13C12 



C13S12 



si3e 



U 



-C23S12 - Ci 2 Si 3 S 2 3e lS C12C23 - Si 2 Sl3S23e I<5 C13S23 
S12S23 - Ci2Si 3 C23e i ' 5 -C12S23 - C23Sl3Sl2e 115 C13C23 

and the neutrino oscillation data at 3cr level [HBj and assuming hierarchical neutrino 



(40) 



masses, 



ATO21 


= (7.09-8.19) x 


10" 5 eV 2 




= (2.18-2.73) x 


10" 3 eV 2 


sin 2 6»i2 


= 0.27- 0.36, 




sin 2 6>23 


= 0.39-0.64, 




sin 2 6> 13 


= 0.092 ±0.06. 





(41) 

We take the leptonic Dirac phase 5 in the U matrix to be zero for which the predicted CP- violation from unitarity 
vanishes irrespective of the values of #13. We have also checked that inclusion of larger values of #13 ~ 8° — 9° [60| do 
not alter our results significantly. Similar results are obtained with 5 ~ tt. 

Taking the light neutrino mass eigen values mi = 0.001 eV, rri2 = 0.0088 eV, — 0.049 eV, and the constructed 
U matrix, we utilize the rj matrix of eg. (1341) for the degenerate case and eq. (|12[) to obtain the nonunitary matrix 
Af. Using eq. (l39|) we also get the matrix. Once the matrices 77 and U are determined as discussed above and in 
Sec. 4, the CP-violating parameters are computed using eq. p^|) . Even though U has no imaginary part because of 
assumed vanishing value of its Dirac phase, CP-violation would arise from the imaginary parts of the corresponding 
components of rj matrix. We also estimate branching ratios for different LFV decay modes using eq. (j!7p . 

For the degenerate case with m,R = 1.1366 TeV we get 
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m Rl2 
(GeV) 


m,R 3 
(GeV) 


Mass eigenvalues /is { 
(MeV) 


50 


3324.7 


(2.4583, 3.23 x 10" 3 , 1.18 x 10~ 6 ) 


100 


1286.5 


(8.0423, 2.60 x 10" 3 , 1.07 x 10~ 6 ) 


500 


1140.7 


(199.37, 5.29 x 10~ 2 , 1.05 x 10~ 6 ) 


1136.6 


1136.6 


(1030.0, 2.72 x 10 _1 , 1.04 x 10 -6 ) 



TABLE V: Mass eigen values of fig signifying masses of singlet fermions predicted by the inverse seesaw in SO(10) 



Ms 



1 0.9932 - 0.0124i -0.1908 + 0.0022i 0.0066 - 0.0033i \ 
-0.1908 + 0.0022i 0.0370- 0.0004i -0.0013 + 0.0006i GeV, 
^ 0.0066 - 0.0033i -0.0013 + 0.0006i 0.00003 - 0.00004i / 



AJ^ = -1.3082 x 10" 6 , 

A ^ 3 = -1-5573 x 10~ 6 , 

Ajf T = 1.5574 x 10~ 6 , 

AJ^ = 1.5572 x 10~ 6 , 

A J T 12 = 4.0144 x 10~ 6 , (42) 



and the branching ratios 



BR(n -> ej) = 2.0025 x 10~ 16 , 
BR(t -> e 7 ) = 2.1586 x 10~ 14 , 

BR{t /i 7 ) = 3.0290 x 10~ 12 . (43) 

Thus we find that in this non-SUSY S'O(IO) model for the degenerate RH neutrino masses, like the SUSY SO(10) 
prediction [29^ . although all the five CP violating parameters are just one order smaller than the corresponding 
parameter in the quark sector where Jckm — (3-051q'2o) x 10~ 5 , there are certain quantitative differences. The 
magnitudes of predicted CP-violations for all the five parameters in the non-SUSY 50(10) model are reduced by 
nearly 50% compared to their corresponding SUSY 50(10) values. 

When compared with the predicted values in SUSY 50(10) [2§] the present results on branching ratios satisfy, 

BR(fi -> e 7 ) susv 3 



BR(p -» 


^)non- susy 


~ 2' 


BR{t 


' /susy 


~ 5, 


BR(t -» 


^^ / )non— susy 


BR{t 


-> W) susy 


2 



BR(t -> M7)non- S usv 3 ' 



(44) 



which can be tested by next generation experiments on LFV decays. 

Our predictions for the partially degenerate RH neutrinos on different elements r\ a $ and their phases are given in 
Table IIIII and those for CP-violating parameters AJ^ and branching ratios are summarized in Table IIVI 

Compared to the predictions in the degenerate case, |?7 M7 -| ~ 10 -4 , S^ T ~ 10~ 6 , for the partially degenerate case 
we find that while |t7 mt | is of the same order, but S^ T ~ 10 -2 , 10 -3 and 10~ 4 for m Rl 2 — 50 GeV, 100 GeV, and 500 
GeV, respectively. These parameters enter into the neutrino oscillation probability in the "golden channel" [H^ , 

P, T ~ 4|^ T | 2 +4 S 2 3 C 2 3 sin 2 (^§^ 

4E 



'Am 2 L N 

4|t7 pt | sin 5^ sin26> 23 sin ( 34 ) . (45) 
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leading to the CP-asymmetry, 



A 



CP 



P — P— 

-A\r] tlT \smS f j.r 

■ nn ■ ( Ami. L 

sin 26*23 sin I ^ 



(46) 



when the hrst term in eq.(|45j) is much smaller compared to the other two terms. Our results in the partial degenerate 
case satisfiy the condition that gives eq. (|4"6l from eq. (|45p . The nonunitarity CP violating effects are predicted to be 
much more pronounced by noting that the strength of the third term in eq. (|45]) is enhanced by 100 — 10, 000 times 
compared to the prediction in the degenerate case. Crucial to this prediction is our constraint eg. (1381) between RH 
neutrino masses which plays an important role in estimating the phase of rj^ T in the partially degenerate case that 
takes into account the increasing behavior of mn 3 for decreasing values of = toa 2 . 

Among other significant differences in the model predictions are Br(/i — > ej) values higher by two orders or by 
one order for = mn 2 = 50 GeV or 100 — 500 GeV while Br(r — > 67) is predicted to be one order lower for 
the RH neutrino masses n is, = to# 2 = 100 — 1180 GeV. Presently the experimental limits on branching ratios are 
Br(/i ->■ ej)< 2.4 x 10" 12 [gfj], Br(r ->■ ej)< 1.2 x 10" 7 and Br(r ->■ (ij)< 4.5 x 10" 8 The projected reach 
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of future sensitivities are up to Br(r — > ej) ~ 10 -9 , Br(r — > (J,-y)~ 10~ 9 , but Br(/i —> 67)^ 10" 

In Table fVl we show predictions of mass eigen values of the /is matrix that signifies masses of three fermion singlets 
Si {i = 1,2,3) for degenerate and partially degenerate cases of RH neutrino masses. These mass eigen values are 
noted to vary starting from the lightest ~ 1 eV to the heaviest ~ 1 GeV which may have interesting phenomenological 
consequences that needs further investigation. It is to be noted that the smallest mass eigen value is also predicted 
directly by the inverse seesaw formula from the TeV scale value of (Md) 33 ~ 100 GeV in a manner similar to the 
Type-I seesaw case. 

We have also examined the consequences of quasidegenerate light neutrino masses expected to manifest through 
tritium beta decay or neutrinoless double beta decay searches. For example with mi = 0.09923 eV, 7712 = 0.09965 
eV, and = 0.1110 eV, which are consistent with neutrino oscillation data, the three eigen values of the resulting 
[is matrix are fig = (30.110 GeV, 1.2 MeV, 20.6 eV) with three pairs of heavy pseudo Dirac neutrinos having almost 
degenerate masses (1151.7, 1121.6) GeV, (1139.5, 1139.5) GeV, and (1136.5, 1136.5) GeV. The predictions for LFV 
decays, CP-violating parameters and the nonunitarity effects are similar to the case of the degenerate pseudo Dirac 
neutrinos with hierarchical light neutrino masses as discussed above. However, the heaviest eigen value of the fermion 
singlet mass matrix increases to fig ~ 30 GeV compared to the corresponding value of /Jjp ~ 1 GeV in the hierarchical 
case of light neutrinos as shown in Table IVl 

The introduction of three additional fermion singlets under 5O(10) needed for the implementation of inverse seesaw 
mechanism may be arg ued to be a limitation of the related GUT models. For that matter, the other 5O(10) models 
ofrefs. [II [3 [51 00 El have utilized these singlets to obtain different interesting results. More recently, the 
superpartners of two out of these three fermion singlets have been demonstrated to be acting as components of inelastic 
dark matter 72j. There is another 5O(10)— based radiative inverse seesaw model which has been designed to explain 
the smallness of the us parameter with the symmetry breaking chain 5O(10) — > SU(5) x U(l) x — > SM x U{\) x 
where more nonstandard fermions and singlets have been found to be necessary [33]. These indicate the popularity 
of 5O(10)-singlet fermion models inspite of the stated limitation. 

In this respect the E$ [36| or SU(3) 3 37 -39| type GUT models do not have this limitation as they contain the 
necessary fermion singlets within their fundamental representations but they also contain a number of additional 
nonstandard fermions. In the absence of any experimental data on the masses of these additional fermions at low 
energies, the determination of the Dirac neutrino mass matrix from fermion mass fits at the GUT scale using the 
bottom-up approach adopted here is not possible. As one major objective of the present work is the prediction on the 



lifetimes of gauge boson mediated proton decay p — > 



on which dedicated search experiments are ongoing [4 1 



This model 



5J7(3) 3 type of GUTs do not serve this objective as the corresponding decays are suppressed [37, 
may be important if, ultimately, proton decay search experiments observe a very large lower limit on the lifetime. 

It has been also argued that because of large size of Higgs representations such as 210^ and 126 h needed in 
5O(10) models employing type-I and type-II seesaw mechanisms, GUT-threshold corrections may give rise to larger 
uncertainties in sin 2 6V predictions and associated mass scale(s) 70] . Counter examples of this result in SO(10) 
having Pati-Salam intermediate symmetry (SU(2)l x SU(2)r x SU(i)c X D = G22AD) with unbroken D-Parity 
have been derived with exactly vanishing GUT-threshold corrections on sin 2 9w as well as on the intermediate scale 
7l|. It has been also shown how threshold corrections can be reduced substantially in other 5O(10) models with 
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naturally plausible constraint that all superheavy components of a SO(10) Higgs representation are degenerate in 



masses [18|,[19|. Noting that the Higgs representation 126// is needed for the implementation of the type-I and type-II 
seesaw mechanisms, and the inverse seesaw needs comparatively much smaller Higgs representation like 16//, the 
possibilities of threshold uncertainties are expected to be correspondingly reduced in our models. In particular our 
minimal Model-I contains neither of the larger Higgs representations 210// and 126//; it requires only the smaller 
representations 45//, 16// and 1Qh ± , 10// 2 . Also in case of the Model-II and its extension, the GUT-threshold effects 
due to superheavy components of Higgs representations 210//, 16// and 10/^ , 10// 2 are expected to be substantially 
reduced compared to the 50(10) model of ref. 19] with G2213D intermediate symmetry because of the absence of the 
large representation 126//. The maximal value of proton lifetime is found to increase by a factor 2(4) due to GUT 
threshold effects in our Model-I (Model-II) over the two-loop predictions. 

Regarding other possibilities of inverse seesaw motivated non-SUSY 5O(10), we find that the minimal single-step 
breaking scenario to the TeV scale gauge symmetry , 5O(10) — > G2113, is ruled out by renormalization group and 
coupling unification constraints. One of the two-step breaking chains, 5O(10) — > G224D —> O2113 gives a low value of 
the unification scale Mjj — 10 14 7 GeV whereas 5O(10) — > G214 — > G2113 also yields almost similar value, My — 10 148 
GeV where we have used SU(2) L x U(1) R x SU(A) C = G 2 u- The third remaining chain, 5O(10) -> G224 -> G2113, 
where D-parity is broken at the GUT scale, gives Mjj = 10 1515 GeV. Thus all the three minimal chains at two-loop 
level are ruled out by the existing lower bouind on proton lifetime (49| . As the large representation 126// is absent 
in these models, the GUT-threshold effects are smaller in the corresponding minimal models than the required 
values to make them compatible with the lower limit on proton lifetime unless the splitting among the superheavy 
components is too large. In view of these, the minimal Model-I turns out to be the best among all possible single and 
two-step breaking minimal models of 5O(10) with the TeV scale G2113 gauge symmetry. 

One of the appealing features which have been noted @ in 5O(10) breaking chains under the category of Model-I 
is that they do not have the cosmological domain wall problem [73j because of spontaneous breaking of D-Parity 
along with the gauge symmetry at the GUT scale. When this criteria is included while searching for equally good 
models, there are only two possible chains with three step breakings and only one chain with four step breaking 
to the TeV-scale symmetry G2113. However, if utilization of large Higgs representations is excluded, the minimal 
Model-I emerges to be unique from among all possible 50(10) breaking chains. Investigation of prospects for these 
longer symmetry breaking chains along with others which is beyond the scope of the present work will be addressed 
elsewhere. 



VIII. SUMMARY AND CONCLUSION 



We have investigated the prospects of inducting TeV-scale inverse seesaw mechanism for neutrino masses into 
nonsupersymmetric 50(10) grand unification and found that it can be successfully implemented with a low-mass 
Z' gauge boson accessible to experimental detection at LHC and planned accelerators. By setting up RGEs in the 
presence of G2213 and G2113 gauge symmetries we have extrapolated fermion masses and mixings to the GUT scale 
using bottom-up approach and determined the Dirac neutrino mass matrix from a successful fit at the GUT scale. 
We have found a relation between the RH neutrino masses which, in the partially degenerate case, predicts the third 
generation RH neutrino mass to increase substantially with the decrease of first or second generation RH neutrino 
masses. Although the predicted branching ratios in the case of degenerate RH neutrinos show less than one order 
variations from the corresponding SUSY SO(10) predictions, in the partially degenerate case, the branching ratio 
Br(/i — ¥ &y) is predicted to be larger by 1 — 2 orders while Br (r — > ej) is predicted to be lower by one order for all 
values of allowed RH neutrino masses. For the nonunitarity matrix element ?7 Mr an important model prediction is its 
enhanced phase 5^ T larger by 2 — 4 orders which is expected to play a dominant role in the experimental detection 
of the nonunitary CP-violation effects at neutrino factories. We have also shown that the models accommodate 
quasidegenerate light neutrino masses relevant for neutrinoless double beta decay or the tritium beta decay searches 
with predictions on the LFV, CP- violation, and nonunitarity effects similar to the case of hierarchical light neutrinos 
and degenerate pseudo Dirac neutrinos while the heaviest mass of the fermion singlets increases from ~ 1 GeV 
to /4 1} ~ 30 GeV. 

Interestingly, the two- loop prediction on proton lifetime in the minimal model (Model-I) turns out to be [r p (p — > 
e + 7r°)] max . = 7 x 10 34 Yrs. which increases by a factor of 2 when GUT threshold effects are included. While providing 
a possibility of verification of the underlying GUT hypothesis, this offers another opportunity for testing the minimal 
model by ongoing search experiments regarding its validity or falsifiability. We have also identified this model to 
be the best among all involving single or two-step breakings of 50(10) to the TeV scale gauge symmetry G2113 
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MODEL 


Symmetry 


04 


hj 
(GeV) 


1, 1 1 


G213 


(-19/6,41/10,-7) 




'199/50, 27/10,44/5 

9/10,35/6, 12 
v 11/10,9/2,-26 , 


) 


I, II 


G2113 


(-3,53/12,33/8, -7) 




/ 8,1,3/2,12 ^ 
3, 17/4, 15/8, 12 
9/2, 15/8, 65/16,4 
\9/2,3/2,l/2,-26/ 




I 


G2213 


(-8/3,-13/6,17/4,-7) 




/ 37/3,6,3/2,12 \ 
6,143/6, 9/4, 12 
9/2,27/4, 37/8,4 

\9/2, 9/2, 1/2, -26/ 




II 


G2213D 


(-13/6,-13/6,17/4,-7) 




/ 143/6,6,9/4,12 ^ 
6,143/6, 9/4, 12 
27/4, 27/4, 23/4,4 
\9/2,9/2,l/2,-26/ 





TABLE VI: One-loop and two-loop beta function coefficients for gauge coupling evolutions in Model-I and Model-II described in the text 

taking the second Higgs doublet mass at 1 TeV 



which is essential for low mass Z' and prominent nonunitarity effects. But if large Higgs representations are excluded 
from symmetry breakings, the minimal Model-I turns out to be a unique model among all pssible 5*0(10) symmetry 
breaking chains. 

Too fast proton decay in another model ( Model-II) has been shown to be evaded by a simple extension where some 
of the predictions on r p should be within the reach of future experiments. On the other hand, if the actual proton 
lifetime is too large, this is also shown to be accommodated in model extensions along with associated nonunitarity 
and lepton flavor violation effects with the prospect of detection of a color octet scalar at accelerator energies. 

In conclusion we find that induction of TeV-scale inverse seesaw mechanism into nonsupersymmetric 5*0(10) pre- 
dicts pronounced nonunitarity and CP-violating effects measurable at accelerator energies and neutrino factories for 
hierarchical as well as partially degenerate spectra of light neutrino masses. In the TeV scale inverse seesaw mechanism 
motivated GUT model, these effects are mainly due to predominance of the Dirac neutrino mass matrix in 5O(10) 
because of its underlying quark-lepton symmetry and this holds even if only an experimentally verifiable low-mass Z 1 
gauge boson is present as one of the smoking gun signatures of asymptotic parity restoration. 
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APPENDIX A 



In the standard notation of two-loop evolution equations for gauge couplings, 



(47) 



the one- and two- loop beta function coefficients are given in Table fVTl We have noted a small contribution of U(1)r x 
U{\)b-l mixing effect [46| especially in the case of Model-I. 
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APPENDIX B 



Each of the two 50(10) models we have considered for inverse seesaw has two types of nonstandard gauge sym- 
metries, G2213 or G2213D and G2113. Here we derive RGEs for running Yukawa and fermion mass matrices from 
which following the earlier approach (fijjjl, we derive RGEs for the mass eigenvalues and mixing angles. We define the 
rescaled /?- functions 



With G2113 symmetry the scalar field $d(2, 1/2,0, 1) through its VEV Vd gives masses to down quarks and charged 
leptons while $ u (2, — 1/2, 0,1) through its VEV v u gives Dirac masses to up quarks and neutrinos. These fields are 
embedded into separate bi-doublets in the presence of G2213 and their vacuum structure has been specified in Sec. 4. 
We have derived the beta functions for RG evolution of Yukawa matrices (Yi), fermion mass matrices (Mi), and the 
vacuum expectation values (v u ,d). The rescaled beta functions are given below in both cases. 



167T 2 /i— — = /3p\. 



(48) 



G2113 Symmetry: 











(49) 



where the beta-functions for VEVs are 



'v. 





(50) 



with 



T u = Tr(3r u ty u + Y }Y U ), T d = Ti(3Y^Y d + Y^Y e ). 



(51) 



The parameters occurring in these equations, and also in eq. (|56[) and eq.f I57p given below are 



a 



3 u 1 ' W n 

2' 2' ' 



C'l 



(9/4,3/4,1/4,8), C\ = (9/4,3/4,9/4,0), 

(0, 0, 1/4, 8), C\ = (0, 0, 9/4, 0), = (9/4, 3/4, 0, 0), (i = 2L, 1R, BL, 3G). 



(52) 
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G 2 2i3 Symmetry: 

Following definitions of Sec. 4 in the presence of left-right symmetry, the rescaled beta functions for RGEs of the 
Yukawa and fermion mass matrices are 

0Y n = (Y u Y^+Y d Yj)Y u + Y u (Y^Y u + Y^Y d )+T u Y u + f 1 Y d -J2c!glY u , 

i 

fty d = (Y d Yj + Y u Y^Y d +Y d (YjY d + Y^Y u )+T d Y d + f 2 Y u -J2CfgfY d , 

i 

Py„ = {Y v Yj + Y e Y%)Y v + Y v {YlY v + Y^Y e ) + T U Y V + f x Y e - £ C l l9 fY„, 

i 

p Ye = (Y e Yl +Y„Y})Y e + Y e {YlY e + Y,}Y„)+T d Y e + f 2 Y„-J2C l ig ?Y ei 



(F u yJ + Y d Yj)M u + M u (YjY u + YjY d ) - £ C?g?M u + Ti tan/3M d 



ftu d = (Y d Yj +Y u Y^)M d + M d (YjY d + Y^Y u ) -J2c!gf}M d + ^-M u , 

i 

ftu D = {Y v uYj + ),), :.\l„ + M D (YJY„ + Y}Y e ) - ^ C l i9 ^M D + T\ tan ftM e 







M r . 



(Y e Yj + Y v Y})M e + M e (Y^Y e + Y}Y U ) - ]T C\g\M e 



tan/3 



-M 



D ■ 



(53) 



where the rescaled beta functions for VEVs ft Vu , (3 Vd are the same as in eq. (|50[) with different coefficients C™ defined 
below and functions T u and T d are the same as in eq. (f!)l"j) . Other two traces entering in this case are 



I\ = Tr(3Y}Y u +YjY v ), 
f 2 = Tv(3Y^Y d + YjY e ). 

The parameters occurring in these equations and also in ea. (|56l) and eg. (1571) given below are 

a = b = 2, a' = b' = 1, 
C? = (9/4,9/4,1/4,8), C\ = (9/4, 9/4, 9/4,0), C? = (0,0, 1/4, 8), 
C\ = (0,0,9/4,0), Cf = (9/4,9/4,0,0), (i = 2L, 2R, BL, 3C). 



(54) 



(55) 



Then following the procedure described in [63j . and using the definition of parameters in the two different mass 
ranges, given above we obtain RGEs for mass eigenvalues and elements of CKM mixing matrix V a p which can be 
expressed in the generalized form for both cases, 

Mass Eigenvalues: 



ft,, 



ftr, 



J2cl q) 9l + ay* + 2b £ \ Vuj ? y * + a 'T±^l £ | Kj f TO . 



j—d,s,b 



j=d,s,b 



To 



tan hi H 51 

A: j—u^c,t j=u,c,t 



to, , i = u, c, £ 



to^ i = d,s,b 



To 



k j=N u N 2 ,N 3 ' j=N 1 ,N 2 ,N 3 



E 



m,-, i = e,(i,T 



S^M 1 ) 2 , 2 , ,7i tan/3 ^ 



TO , 



to,, i = N u N2,N s 



(56) 
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CKM Matrix Elements: 



7— ii,c,t;7^a 



, Titan/3 y , 



b mi 



mz 



V a 



1 " 

r y=d,s,b; r Y^/3 



r 2 



tan/ 



(V^M u V) lP 



b mz) + mj 



(57) 



Then using third generation dominance, the beta functions for all the 9 elements are easily obtained for respective 
mass ranges where in addition to the parameters in the respective cases in ea. (l52[) and eq. (|55[) . a' = b' = in the 
mass range M R o — ► M R + with G2113 symmetry, but a' = b' = 1 in the mass range M R + — > Mjj with G2213 or G2213D 
symmetry and, in the latter case, the nonvanishing traces T\ 2 are easily evaluated in the mass basis. 
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